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ABSTRACT 

A mathematical  model  which  was  developed  by  Monaghan  at  DRES 
was  extended  to  predict  the  penetration  of  vapour  through  clothed  skin 
for  an  initial  liquid  load  on  or  in  the  clothing.  The  model  and  its 
associated  computer  program  along  with  some  sample  calculations  are 
described  tn  this  report. 
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decomposition  rate  (min-1) 

diffusive  and  absorptive  coefficients  for  slice  i 

saturated  vapour  concentration  In  air  (mg  cm-3) 
computer  language  for  At 

general  boundary  term  introduced  to  account  for  vapour  evolution 
In  the  inward  direction  from  a vapour-liquid  interface  (mg  cm"2) 

same  as  g,  except  for  vapour  evolution  in  the  outward  direction 
(mg  cm"2) 

subscript  for  time  increment  number 
number  of  slices  in  a clothing  layer 
mass  per  unit  area  in  slice  i (mg  cm-2) 

total  amount  of  surface  liquid  lost  to  vapour  evolution  (mg  cm"2) 
diffusive  resistance  of  air  (min  cm"1) 
diffusive  resistance  of  slice  1 (min  cm"1) 
time  (min) 

time  increment  size  (min) 
total  vapour  loss  to  air  (mg  cm"2) 

absorptive  capacity  of  clothing  or  skin  per  unit  area  of  slice  1 
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1.  INTRODUCTION 


A mathematical  model  (Monaghan)  was  developed  at  DRES  to  predict 
the  penetration  of  liquid  or  vapour  through  clothed  skin.  The  essential 
features  of  the  model  are  described  as  follows: 

The  problem  of  agent  penetration  through  clothed  or 
bare  skin  is  treated  as  a case  of  one-dimensional  diffusion  in 
a multilayer  system  where  absorptive  and  diffusive  properties 
may  vary  from  one  layer  to  another  and  agent  may  decay  in  the 
system  as  a result  of  decomposition.  The  flow  of  agent  in  such 
a system  is  analogous  to  the  flow  of  heat  through  a series  of 
slabs  or  the  flow  of  current  in  a resistance-capacitance  net- 
work when  the  resistance  and  capacitors  are  very  numerous 
(Pattle  and  Monaghan,  1964). 

Experimental  evidence  suggests  that  the  simplest  model 
for  the  skin  is  a system  in  which  two  absorbent  layers  overlay 
a sink  for  agent,  and  only  the  layer  next  to  the  sink  has  a 
significant  diffusional  resistance.  The  two  layers  are  broadly 
identifiable  with  the  homy  surface  layer  and  the  transitional 
layer  of  the  epidermis  and  are  so  named  in  the  model.  The  live 
epidermis  and  the  dermis  are  ignored  since  both  present  no 
appreciable  barrier  to  agent  penetration  and  the  latter  is  per- 
fused by  the  blood  which  acts  as  the  sink  for  agent.  Similarly 
a clothing  layer  can  be  represented  simply  by  a layer  with  a 
significant  absorbence  and  diffusional  resistance.  Decom- 
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position  of  agent  has  been  observed  in  both  skin  and  cloth 
and  is  approximated  by  a first  order  reaction. 

An  electrical  analogue  was  chosen  as  the  basis  of 
the  model.  The  following  table  gives  the  analogues  of  elec- 
trical charge,  potential,  capacitance  and  resistance  used  in 
the  model  for  a clothing  or  skin  layer  of  thickness  l,  with  a 
diffusion  coefficient  D,  that  has  absorbed  a mass  of  agent 
per  unit  volume  C. 


ANALOGUE 

SYMBOL 

UNITS 

CHARGE 

Mass  absorbed  per  unit 
area  = Cl 

m 

mg/cm2 

POTENTIAL 

Equilibrium  concentration, 
vapour  in  air,  Cg  = C/B 

Ce 

mg/cm3 

CAPACITANCE 

Absorptive  capacity  per 
unit  area  =3 l - m/Cg 

V 

cm 

RESISTANCE 

Diffusional  Resistance  = 

1/  DB  = l2/0\l 

R 

min/cm 

The 

coefficient  B and  the  capacitance  V 

are  assumed 

constant  for  a given  agent  and  layer,  but  may  vary  from  layer 
to  layer.  This  assumption  implies  that  V is  independent  of 
concentration. 


The  problem  of  agent  penetration  through  skin  is 
solved  by  deriving  equations  of  flow  from  an  analogous 
electric  R-C  network,  in  which  each  layer  is  divided  into 
a sufficient  number  of  uniform  slices  normal  to  the  flow 
of  agent  to  provide  the  required  accuracy  of  solution. 

The  equations  of  flow  form  a set  of  first  order  differential 
equations  for  each  case  considered  and  are  solved  numerically 
by  means  of  a digital  computer. 

The  computer  program  (McPherson)  was  written  to  calculate  an 
approximate  solution  to  the  differential  equations  and  to  set  up  the 
initial  and  boundary  conditions  of  the  simulations. 

In  this  program  the  model  is  subdivided  into  five  submodels, 
two  of  which  simulate  cases  where: 

a)  liquid  is  deposited  on  the  top  of  the  outside  layer  of 
clothing  on  a clothed  person;  and 

b)  liquid  is  pressed  through  the  clothing. 

In  the  first  of  these  submodels,  the  initial  liquid  loading  is  specified,  and 
in  the  second,  the  initial  liquid  loading  in  the  cloth  is  specified,  but 
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both  loadings  cannot  be  specified  at  once.  Also,  the  second  submodel  and 
its  corresponding  computer  program  did  not  produce  reasonable  results. 

The  objective  of  this  paper  is  to  describe  an  improved  model 
which  can  handle  the  generalized  initial  condition  of  liquid  on  and  in 
cloth. 

2 . FEATURES  OF  THE  MODEL 

The  penetration  model  was  developed  in  order  to  better  understand 
the  hazards  to  people  posed  by  toxic  chemicals  in  the  liquid  phase.  In 
other  words,  this  study  deals  with  an  attempt  to  describe  mathematically 
the  physical  processes  occurring,  with  time,  when  liquid  droplets  impact 
on  the  clothing  of  a person,  or  when  liquid  is  forced  into  the  clothing  by 
applying  pressure,  e.g.  by  sitting.  The  passage  of  the  contaminant  through 
clothing  and  skin  to  the  blood  stream  is  a time  dependent  diffusion  process. 

In  addition  to  the  above,  two  other  provisions  are  part  of  the 
model.  If  the  challenge  is  to  a clothing  ensemble,  then  provision  is  made 
for  the  removal  of  this  contaminated  clothing  at  some  point  in  time.  This 
is  accomplished  by  zeroing  the  mass  that  exists  in  each  of  the  clothing 
layers  at  that  time.  Similarly,  the  model  will  simulate  the  decontamination 
of  the  outer  clothing  surface  by  zeroing  the  surface  contaminant. 

In  its  present  form,  the  model  can  handle  up  to  four  layers  of 
clothing  over  skin.  Each  layer  of  clothing  is  assumed  to  be  homogeneous. 

The  skin  is  divided  into  three  distinct  regimes,  the  horny  surface  layer, 
the  transitional  layer  of  the  epidermis,  and  the  blood  stream  (a  sink). 

3.  CALCULATION  OF  VAPOUR  PENETRATION 

a)  Basic  Mathematical  Relationships 

Penetration  is  described  by  a one-dimensional  diffusion  equation 
in  a multi-layer  system,  i.e.  clothing  and  skin,  where  absorptive  and  dif- 
fusion properties  may  vary  from  one  layer  to  another  and  mass  absorbed  may 
decay  with  time  in  the  system  as  a result  of  decomposition. 

Each  layer  of  the  system  is  divided  into  several  slices,  the 
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number  of  which  is  chosen  to  provide  adequate  computing  accuracy. 

A general  expression  for  flow  of  agent  vapour  into  a slice  can 
be  obtained  by  applying  Kirchoff's  law  at  a node  in  the  analogous  re- 
sistance capacitance  network  (Fig.  1).  Note  that  Kirchoff's  law  states 
that  at  any  junction  in  the  circuit  the  total  current  flowing  toward  the 
junction  must  equal  the  total  current  flowing  away.  Considering  flow  into 
slice  i we  obtain: 


c*mi  1 

. M 

i 2 J 

S+l  l"i\  2 

dt 

lVi 

vi ) 

,Ri-l+Ri  ' 

\Vl  Vi/Ri+Ri+1 

where 

rru  is  the  mass  per  unit  area  in  slice  i (mg  cm"2) 

R.j  is  the  diffusive  resistance  of  slice  i (min  cm-1) 
V.j  is  the  absorptive  capacity  per  unit  area  of  slice 
a is  the  decomposition  rate  (min-1) 


(Eq.  1) 


(cm) 


A simple  implicit  difference  method  was  used  to  solve  the  dif- 
fusion equation  (Eq.  1).  If  we  denote  the  mass  per  unit  area  of  skin  slice 

i at  time  jAt  by  m.  ■,  and  m.  . , at  the  next  increment  of  time,  then  Eq.  1 
* > J * * J+  * 

can  be  written  in  the  following  discretized  form: 


mi,j 


m 


i , j+1 


fmi-l,j+l 

1 2At 

fmi+l,j+l 

m.  ..,N 
i.J+l 

\ 2At 

l v, 

' »,  i 

'Ri-1+Ri 

[ v, 

' V,  , 

' Ri+Ri+1 

Now  let  Ai  = 


-2At 


VRi  + ,w 


Bi  ■ 1 + 


ci 


’l_ 

2 _ . + _ 2 

i.  a 

v. 

[Ri.l  + Ri  Ri  + Ri+ly 

• a 

l 

_ 

At 


-2At 


Vi(Ri-l  +Ri) 


(Eq.  2) 
(Eq.  3) 

(Eq.  4) 

(Eq.  5) 


Then  Eq.  2 can  be  written: 

* chi"i+i,j+i  * mi ,j  6> 

Now  the  general  expression  for  vapour  diffusion  in  a skin-clothing  system 
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with  n slices  is  s‘ven  in  the  following  matrix  form 


B,  C2  0 0 

A1  B2  C3  0 

o a2  b3  c4  0 .... 

_ 

m,  . , . 

1 ,J+1 

mi  J 

A B , C 

n n-1  n 

0 ....  0 A . B 

n-1  n_ 

n,j+l  J 

_ ""’J  _ 

(Eq.  7) 


b)  Boundary  and  Initial  Conditions 

While  the  basic  mathematics  of  the  model  appears  quite  simple, 
the  elegance  of  it  is  attained  through  the  control  of  the  various  boundary 
and  initial  conditions. 

The  initial  conditions  consist  of  the  liquid  loading  on  the 
outer  clothing  surface,  and  the  magnitude  and  distribution  of  the  liquid 
loading  within  each  layer  of  clothing.  The  model  does  not  allow  gaps  in 
the  liquid  loading.  All  the  liquid  loads  must  be  contiguous. 

Sperial  conditions  occur  at  the  inside  and  outside  boundaries  of 
the  clothed  skin  system  and  at  liquid-vapour  interfaces  (Fig.  2).  These 
are  accounted  for  by  modifying  the  corresponding  equations  in  the  general 
expression  (Eq.  7)  as  follows: 

(1 ) First  equation  in  the  system 

The  first  slice  in  the  clothed  skin  system  represents  the  sys- 
temic sink,  the  potential  and  resistance  of  which  are  zero. 

Applying  Kirchoff's  law  for  flow  into  the  slice  gives: 


dm-|  2m2 
dt  V2R2 

Discretizing  Eq.  8 in  the  same  way  as  for  Eq.  1 gives: 


(Eq.  8) 


ml,j+l  + C2m2, j+1  = ml,j  (Eq’  9) 
Comparison  to  the  general  expression  given  by  Eq.  4 and  Eq.  6 
then  gives: 
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Blml,j+.  + C2m2, j+1  = mij  + gF  (Eq*  10) 

where  B-|  = 1 

gF  is  a general  boundary  term  which  was  introduced  to  account 
for  vapour  evolution  in  the  outward  direction  from  a liquid- 
vapour  interface.  In  Eq.  10,  gF  = 0. 

(2)  Inside  of  skin 


The  first  inside  slice  of  skin  is  represented  by  the  second 
equation  in  the  system.  Noting  that  the  resistance  of  the 
first  slice  is  zero  since  it  represents  the  systemic  sink, 
application  of  Kirchoff's  law  for  flow  into  the  second  slice 
gives: 


+ am.. 


^2+B3 


(Eq.  11) 


Then  discretizing  as  for  Eq.  1 gives: 

B2m2,j+1  + C3m3,j+1  = m2,j  (Eq.  12) 

Comparison  to  the  general  expression  given  by  Eq.  6 then  shows 
A1  = 0. 

( 3 ) Inside  liquid-vapour  interface 

The  last  inside  slice  without  liquid  is  designated  i = nk 
(Fig.  2).  The  potential  for  flow  from  the  adjacent  liquid  con- 
taining slice  is  given  by  the  saturated  vapour  concentration, 

C . Application  of  Kirchoff's  law  for  flow  into  slice  nk  gives 


d"nk  _ j 

1 "nk-l 

.%k) 

2 , L mnk  ) 

i 2 

dt  j 

( Vnk-1 

Vnkj 

Rnk-1+Rnk  \ S Vnk  J 

lRnk 

- amnk  (Eq.  13) 


Discretizing  as  for  Eq.  1 gives: 

flnk-l"nk-l,j*l  * Bnkmnk,j+1  = "nk  * B (Eq.  14) 

where  g is  a general  boundary  term  which  was  introduced  to 
account  for  vapour  evolution  in  the  inward  direction  from  a 
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liquid-vapour  interface.  In  Eq.  14,  g = 


2CsAt 


(Eq.  15) 


Comparison  to  the  general  expression  given  by  Eq.  6 gives: 


B . = 1 + 
nk 


^nk  \Rnk-l+Rnk  Rnk 


+ a At 


(Eq.  16) 


For  the  special  case  of  liquid  in  the  clothing  slice  adjacent 
to  the  skin,  Eq.  15  and  Eq.  16  show  that  g and  Bnk  become  in- 
finite because  the  outer  skin  slice  has  no  resistance.  This 
is  not  a reasonable  physical  representation.  To  adjust  this, 
the  assumption  was  made  that,  over  a finite  time  interval,  the 
clothing  dries  from  the  inside  so  that  there  is  a non-zero  re- 
sistance to  vapour  transport.  This  resistance  was  assumed  to 
increase  from  0 to  Rn(c+^  with  time  as  the  clothing  dries.  Then 
the  average  value  of  the  resistance  was  used  in  the  discreti- 
zation, so  that  for  the  outer  skin  slice  Eq.  15  and  Eq.  16 
become: 


2CsAt 


(Eq.  17) 


and  Bnk  - 1 + 


Vnk  \Rnk-l 


+ a At 


(Eq.  18) 


(4)  Outside  liquid-vapour  interface  within  clothinc 


The  first  outside  slice  without  liquid  is  designated  i = nl 
(Fig.  2).  As  for  the  previously  mentioned  interface,  the  poten- 
tial for  vapour  transport  from  the  adjacent  slice  containing 
liquid  is  the  saturated  vapour  concentration,  C$.  Application 
of  Kirchoff's  law  for  flow  into  slice  at  gives: 


V jj_+  V+i  _ V 


Vn£+1  V/  Rn£+Rn£+1 


- am  , (Eq.  19) 


Discretizing  as  for  Eq.  1 gives: 
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BlV,j+l  + Cn£+lV+l,j+l  = m 


nl,j 


i + 9p 


where  gF  = 


2C$  At 


R 


(Eq.  20) 
(Eq.  21) 


nl 


Comparison  to  the  general  case  by  referring  to  Eq.  4 and  Eq.  6 
shows  that  the  value  B is  given  by: 


Bn£  = 1 + 


" 1 

2 + 2 11 

V 

\Rnl  Rn£+Rn£+1  ) 

At 


(Eq.  22) 


(5)  Outside  clothing  slice 

The  outside  clothing  slice  is  designated  by  i = nn  (Fig.  2). 

The  vapour  transport  equation  for  this  slice  is  the  last  in 
the  system  of  n equations.  Liquid  may  or  may  not  be  present 
on  the  surface.  This  must  be  accounted  for  in  the  equations 
for  this  slice.  Also,  there  is  the  special  case  of  only  one 
slice  without  liquid,  which  must  be  treated  separately. 

(5a)  Two  or  more  slices  without  liquid 

1.  Liquid  on  the  surface.  Application  of  Kirchoff's 
law  for  flow  into  slice  nn  gives: 


dni 

nn 

ini  . hi 
nn- I nn 

2 , 

fc  . Hlnn.1 

2 

dt 

1 V V 

( nn-1  nn  j 

Rnn-1+Rnn 

lS  'n  n| 

Rnn 

- am 


Discretizing  as  for  Eq.  1 gives: 
Ann-lmnn-l,j+l  + Bnnmnn,j+1  = mnr.,j  + 9 


where  g = 


2CsAt 


nn 


nn 

(Eq.  23) 
(Eq.  24) 
(Eq.  25) 


Comparison  to  Eq.  4 and  Eq.  6 gives  the  modified  value  of 
B as  follows: 


B 


nn 


1 + 


■ 

1 

f 2 f 2 1 

J.  a 

V 

Rnn  R 

\ nn-1  nn  nn  ^ 

T a 

1 nn 

At 


(Eq.  26) 
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2.  No  liquid  on  the  surface.  In  Eq.  23  there  is 
zero  potential  in  place  of  C$  and  resistance  for  flow 
rate  into  the  slice  from  the  outside  is  changed  from 
R R 

to  + R where  R is  the  diffusive  resistance  of  air. 

2 2 d ^ 


Then  g = 0 


and  B„  = 1 + 
nn 


(Eq.  27) 


V R ,+R 
nn  nn-l  nn 


+ V a At  (Eq.  28) 

R +2R, 
nn  a/ 


(5b)  Special  case  when  only  one  outside  slice  does  not 
contain  liquid 

Application  of  Kirchoff's  law  for  flow  of  vapour  into 
slice  nn  gives  the  following  equation: 


— C - -£-+  o - 

dt  \$  Vnn/Rnn  l V„r,  \ 

Discretizing  as  for  Eq.  1 gives: 

B m„„  ,xi  = rn  ^ + gF  + g 
nn  nn,j+l  nn,j  3 3 

where  g = 0 


and  gF  = 


2c$At 


(Eq.  29) 

(Eq.  30) 
(Eq.  31) 

(Eq.  32) 


Comparison  to  Eq.  4 and  Eq.  6 gives: 


B = 1 + 
nn 


V nn  Rnn  Rnn+2Ra  j 
nn  \ nn  nn  a / 


(Eq.  33) 


c)  Complete  System  of  Equations  in  Matrix  form 

The  equations  of  vapour  transport  can  now  be  written  in  the  form 
of  Eq.  7 with  all  the  modifications  to  account  for  the  various  boundary 
conditions. 

Considering  boundary  conditions  (1)  to  (3)  the  equations  of 
vapour  transport  in  the  inside  slices  are  written  as  follows: 
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— —1 

““  * — 

— i 

B-j  Cg  0 . , . . 0 

A1  B2  ^3  B 

o a2  b3  c4  0 

ml , j+1 
• • 

• • 

• • 

"l.j 

• • 

• • 

+ 

gE 

0 

• • 

• • 

..  ..  0 A|ik_1  Bn|(_1  Cnk 

• • 

• • 

• • 

• • 

0 

_°  **  **  0 Ank-1  Bnk_ 

L_mnk,j+1_ 

_mnk,i 

g 

where 


A1  = 0 

B1  = 1 
gF  = o 


g = 


2CsAt 


'nk 


Bnk  = 1 + 


' i 

2 , 2 

*■ 

X a 

_Vnk 

(Rnk-1+Rnk  Rnk] 

“ a 

m 

At 


(Eq.  34) 


(Eq.  35) 
(Eq.  36) 
(Eq.  37) 

(Eq.  38) 

(Eq.  39) 


except  when  the  slice  adjacent  to  the  skin  contains  liquid.  Then: 

2C  .At 

g = — (Eq.  40) 
Rnk+1 


and  Bnk  = 1 + 


" 1 

f 2 2 \ 1 

> 

^Rnk-1  Rnk+1  / 

At 


(Eq.  41) 


which  is  identical  to  Eq.  4 because  Rnk  = 0. 

Considering  boundary  equations  (4)  and  (5)  the  matrix  system  for 
vapour  transport  in  the  outside  slices  is  written  as  follows: 


"■  — 

p — 

_ — 

Bn l Cn£+1  0 • * 0 

V Bn*+1  Ca£+2  0 * * 

0 ••  ••  ••  •• 

• • • • ••  ••  • • 

V,j+i 
• • 

• • 

• • 

v,j 
• • 

• • 

• • 

+ 

gF 

0 

• • 

• • 

• • ••  ••  ••  •• 

• • 0 A •»  B i C 

nn-i  nn-i  nn 

• • 

• • 

• • 

• • 

0 

0 • • 0 A -I  B 

_ nn- 1 nn_ 

%in,j+l_ 

v,j_ 

g 

(Eq.  42) 
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where  gF  = 


2C$At 


'n  l 


' 1 

2 + 2 

y 

l Rae  W+l  j 

= 1 + 


1 ) Liquid  on  Surface 


At 


/II 

(Eq.  43) 

(Eq.  44) 


9 = 


2C  At 
s 


nn 


L = 1 + 
nn 


^nn  lRnn-l+Rnn 


+ a 


nn 


At 


(Eq.  45) 

(Eq.  46) 


2)  No  Liquid  on  Surface 
g = 0 


B = 1 + 
nn 


1 


V [R  i+R  R +2R 
nn  I nn- I nn  nn  a 


+ a 


At 


(Eq.  47) 
(Eq.  48) 


When  there  is  only  one  outside  slice  which  does  not  contain 
liquid  the  matrix  system  of  Eq.  42  reduces  to  a single  equation  as  follows: 


j.t  = m„„  . + gF  + g 
nn  nn,j+l  nn,j  3 3 


where  gF  = 


2C$At 


R_r 

nn 

and  g = 0 

B - ! + 

1 

f 2 + 2 U al 

nn 

Vnn  1 

» 

J 

At 


(Eq.  49) 
(Eq.  50) 

(Eq.  51) 
(Eq.  52) 


When  there  is  no  liquid  at  all  in  the  clothing,  the  equations 
for  vapour  transport  can  be  expressed  by  one  matrix  system  as  follows: 
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~B1 

C2 

0 

• • 

• • 

0 

ml,j+l 

’ "l.j 

gF 

A1 

B2 

C3 

0 

• ■ 

• • 

• • 

0 

0 

A2 

B3 

B4 

0 

• • 

• • 

• • 

• • 

• • 

• • 

• • 

• • 

• • 

• • 

= 

• • 

+ 

• • 

• • 

• • 

• • 

• • 

• • 

• • 

• • 

• • 

• • 

A i 

B n 

c 

o 

nn-1 

nn-1 

nn 

0 

• • 

• • 

0 

Ann-1 

1 

(Z 

C 

cc 

%n,j+l 

g 

(Eq.  53) 


where 


D 


2) 


A1  = 0 
B]  = 0 
gF  = 0 

Liquid  on  Surface 


9 = 


2C$At 


nn 


' l 

2 , 2 ' 

4.  a 

V 

tRnn-l+Rnn  Rnn  ) 

▼ Q 

| nn 

B = 1 + 
nn 


No  Liquid  on  Surface 
g = 0 


At 


B = 1 + 
nn 


" 1 

2 2 ^ 

■ 

+ a 

V 

R ,+R_n  Rnn+2R 
k nn-l  nn  nn  a / 

nn 

- 

At 


(Eq.  54) 
(Eq.  55) 
(Eq.  56) 


(Eq.  57) 
(Eq.  58) 


(Eq.  59) 
(Eq.  60) 


All  three  matrices  given  by  Eq.  34,  Eq.  42  and  Eq.  53  are  of 
the  same  font)  and  can  be  written  in  matrix  notation  as  follows: 

A Mj+]  = Mj  + G (Eq.  61) 

where 

Mj  is  the  mass  of  vapour  per  unit  area  in  each  slice  of 

the  system  at  a given  point  in  time  jAt 

Mj+i  is  the  vector  of  mass  per  unit  area  after  the  next 

time  increment 

A is  the  tridiagonal  matrix  of  coefficients  of  the 
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absorptive  and  diffusive  properties  of  the  system 

-► 

G is  the  vector  of  boundary  conditions 


For  a given  set  of  initial  conditions,  the  solution  to  the  pro- 
blem, which  can  readily  be  obtained  by  the  Gauss  elimination  method 
(Westlake),  is  then: 

Mj+i  = A-1  (Mj  + G)  (Eq.  62) 

for  the  two  matrices  given  by  Eq.  34  and  Eq.  42  when  the  clothing  contains 
liquid  and  for  the  matrix  of  Eq.  53  when  the  clothing  contains  only  vapour, 
except  when  only  one  outside  slice  does  not  contain  liquid.  Then  the 
solution  for  this  slice  is: 


_ mnn,.i  + 9F 
mnn,j+l ^ 


B 


nn 


(Eq.  63) 


4.  AUXILIARY  RELATIONSHIPS 
a)  Loss  of  Liquid 

As  vapour  is  transported  away  from  liquid-vapour  Interfaces,  the 
amount  of  liquid  in  the  liquid  containing  slice  must  be  reduced  to  balance 
it.  In  addition,  there  is  a loss  due  to  decomposition,  which  also  occurs 
in  internal  slices  within  the  clothing  layers  containing  liquid.  The 
liquid  initially  on  the  clothing  surface  is  gradually  lost  due  to  evapo- 
ration to  the  air,  and  when  the  slice  of  clothing  next  to  the  surface  con- 
tains no  liquid,  there  is  an  additional  loss  due  to  vapour  transport  into 
the  clothing. 

The  liquid  losses  were  accounted  for  in  each  time  step  along 
with  the  vapour  transport,  previously  described,  as  follows. 

(1 ) First  inside  slice  containing  liquid 

The  first  Inside  slice  containing  liquid  is  designated  nk  + 1 

(Fig.  2).  The  loss  of  liquid  from  this  slice  was  determined 

using  the  electrical  analogy  (Fig.  1)  for  vapour  transport. 

Application  of  Kirchoff's  law  for  flow  into  slice  nk  + 1 gives: 
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‘ aCsVnk+l 


Discretizing  for  time  step  jAt  to  ( j+1 ) At  then  gives: 


(Eq.  64) 


nnk+l , j+1  = mnk+l ,j  + 


2mnk.jAt  2CsAt 


VnkRnk 


aCsVnk+]At  (Eq.  65) 


There  is  a special  case  when  the  first  inside  slice  containing 
liquid  is  adjacent  to  the  skin.  The  value  of  Rnk  is  zero  for 
the  outer  skin  slice,  which  gives  rise  to  terms  of  infinite  value 
in  Eq.  65.  This  was  adjusted  in  the  same  manner  as  for  vapour 
transport,  described  in  section  3. b) ( 3) . Thus,  for  slice  nk+1, 
(Eq.  65)  is  replaced  by  the  following  equation: 


nnk+l , j+1  = "Vik+l.j 


,At  2C$At 


VnkRnk+l 


aEsVnk+lAt  (Eq«  66) 


(2)  First  outside  slice  containing  liquid 

The  first  outside  slice  containing  liquid  is  designated  n£-l 
(Fig.  2).  Application  of  Kirchoff's  law  for  flow  into  this  slice 
gives: 


' aCsVnk+l 


(Eq.  67) 


Discretizing  for  time  step  jAt  to  (j+1  )At  then  gives: 

2m  » .At  2C  At 

V-l,j+l  = V-l,j  + v n ' aCsVe-lAt  (Eq‘  68 } 

W V 

There  is  also  a special  case  for  the  first  outside  slice  containing 
liquid.  This  occurs  at  the  outside  slice  of  clothing,  which  is 
designated  i * nn  (Fig.  2).  When  there  is  no  liquid  on  the  surface 
Kirchoff's  law  for  flow  into  slice  nn  gives: 


- aC  V At 
s nn 


(Eq.  69) 
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Discretizing  gives: 


C At 


itt  i “ m • 
nn,j+l  nn  j j 


■ aCcVnnAt 

s nn 


(Eq.  70) 


When  there  is  liquid  on  the  surface,  Eq.  70  reduces  to: 


nr  . - aC  V At 

nn,j+l  nn,j  s nn 

(3)  Middle  slices 


(Eq.  71) 


The  middle  slices  in  the  clothing  are  designated  i,  where  nk  + 1 
< i < n£  - 1 (Fig.  2).  Since  the  reduction  of  liquid  in  these 
slices  is  by  decomposition  only,  Kirchoff's  law  for  flow  into 
one  of  these  is: 


— = - aC  V.  nk  +1  < i < nl  - 1 (Eq.  72) 

dt  s 1 

Discretizing  gives: 

mi.0+l  + mi , j " aCsVi  nk  + 1 < i < n£  - 1 (Eq.  73) 

Note  that  the  rate  of  loss  by  decomposition  throughout  the 
model  is  assumed  to  be  aC$V.  whenever  m./V..  s Cs.  Also  note 

that  Cs  is  only  used  as  a potential  whenever  m../!^  > C$. 

(4)  Only  one  slice  containing  liquid 

When  there  is  only  one  slice  in  the  clothing  which  contains  liquid 
nk  + 1 = nl  - 1 (Fig.  2).  The  loss  to  the  inside  is  given  by  Eq.  65 
or  Eq.  66,  and  the  loss  to  the  outside  is  given  by  one  of  Eq.  68, 

70  or  71.  However,  in  order  not  to  apply  the  decomposition  twice 
the  quantity  aCs  Vn£+1  At  was  added  to  the  equation  for  the  outside 

slice  whenever  nk  + 1 = n£  - 1. 

(5)  Loss  of  liquid  from  outside  clothing  surface 

When  there  is  no  liquid  in  the  outer  slice  of  clothing,  Kirchoff's 
law  for  the  flow  of  vapour  away  from  the  liquid  on  the  clothing 
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surface  gives: 


(Eq.  74) 


where  Q is  the  total  amount  of  surface  liquid  lost  to  vapour 
evolution. 


Discretizing  gives: 


Qj + 


2Csat 


nn 


2mAt  C.At 
nn  s 

VnnR„n  Ra 

nn  nn  a 


(Eq.  75) 


When  there  Is  liquid  In  the  outer  slice  of  clothing  Eq.  75  reduces 
to: 

C At 

Vi  ■ Qj  * 7~  <E<-  76> 


b)  Total  Loss  of  Vapour  to  Air 


The  Initial  liquid  on  and  In  clothing  Is  lost  to  air,  to  vapour 
In  clothing,  to  the  systemic  sink  and  to  decomposition.  After  the  liquid 
Is  gone  from  the  surface,  there  is  a further  vapour  loss  to  air  from  the 
clothing.  Now,  If  the  vapour  loss  to  air  Is  accounted  for,  it  is  possible 
to  calculate  the  loss  due  to  decomposition,  because  this  Is  the  mass  de- 
ficiency In  the  sums  of  the  vapour  and  liquid  In  the  total  system.  Alter- 
natively, If  the  decomposition  Is  set  to  zero,  the  calculation  of  vapour 
transport  can  be  checked  at  each  time  step  by  observing  that  the  total 
mass  In  the  system  must  be  constant. 


The  loss  of  vapour  to  air  was  accounted  for  in  each  time  step 
as  follows: 

(1 ) No  liquid  in  the  outside  clothing  slice  or  on  the  surface 


Application  of  Klrchoff's  law  using  the  electrical  analogy  (Fig.  1) 
for  vapour  transport  to  air  from  the  outside  clothing  slice  gives: 


dU 

dt 


— 2 \ 
V R+R, 
nn  nn  a / 


(Eq.  77) 
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where  U is  the  total  vapour  loss  to  air. 


Discretizing  gives: 


2m  At 

Ui+1  = U.  = ^ 

' ° W2Ra> 


(Eq.  78) 


(2)  Liquid  in  the  outside  clothing  slice  or  on  the  surface 


Kirchoff's  law  gives:  — = — 

dt  R 


(Eq.  73) 


Discretizing  gives: 


•Vi  ■ + 


(Eq.  74) 


c)  Maximum  Allowable  Size  of  Time  Increment 

The  chosen  time  increment  for  the  calculation  of  vapour  transport 
must  be  such  that  the  mass  in  a slice  never  becomes  negative,  otherwise 
solutions  become  unstable.  Experience  has  shown  that  this  can  occur  for 
liquid  reduction  in  slice  nk+1  as  calculated  from  Eq.  65  when  the  increment 
is  too  large.  Analysis  of  Eq.  65  then  provided  a relationship  for  predeter- 
mining an  approximate  upper  limit  for  the  time  step  size.  This  was  done  in 
the  following  way: 

For  mnk+l  j+1  s ® Eq.  65  we  must  have: 

4t  ( rf-  ■ r aCsVnk+l  s «Wl , j (E"-  81 ) 

\vnknk  Knk  / 

Liquid  must  be  present  whenever  Eq.  65  is  used.  Therefore: 

mnk+l , j * CsVnk+l  (Eq.  82) 

Substitution  of  Eq.  82  into  Eq.  81  and  rearranging  gives: 


CsVnk+l 


(Eq.  83) 


VnkRnk 


A simplified  but  slightly  stronger  restriction  on  the  maximum 


i 
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size  of  the  time  increment  is  obtained  by  setting  the  second  term  in 
the  denominator  of  Eq.  83  to  zero. 

Thus  At  < -k— k-+1-  (Eq.  84) 

2+aR  . V ... 
nk  nk+1 

5.  COMPUTER  PROGRAM 

A computer  program  was  written  to  solve  the  problem  of  vapour 
transport  by  means  of  the  mathematical  model  described  in  the  previous 
sections.  A listing  of  the  program  along  with  one  set  of  output  results 
is  shown  in  Appendix  A.  The  listing  is  annotated  in  detail  to  describe 
the  function  of  each  part  of  the  program  and  the  input  data  required. 

The  program  is  written  in  the  basic  Fortran  IV  language  of  the 
DRES  1130  model  2C  computer  which  has  16K  words  of  core  storage.  The 
program  in  its  present  form  requires  approximately  9K  of  core  storage. 

It  will  allow  up  to  four  layers  of  clothing  to  be  simulated  over  skin 
with  a total  of  fifty  slices  in  the  system.  The  skin  is  usually  divided 
into  seven  slices,  which  leaves  forty  three  slices  for  the  clothing  layers. 

6.  RESULTS 

The  computer  program  was  applied  to  compare  the  systemic  dose 
for  various  distributions  of  a 1.2  mg  cm-2  initial  load  on  and  in  two 
layers  of  clothing.  The  absorptive  and  diffusive  properties  used  were 
chosen  for  test  purposes  only,  and,  in  the  author's  knowledge,  do  not 
represent  any  real  liquid  and  clothed  skin  system.  The  resistance, 
capacitance  and  saturated  vapour  concentration  used  are  those  shown  in 
the  sample  computer  output  (Appendix  A).  The  calculations  were  done  for 
two  different  windspeeds  to  show  the  effect  of  varying  windspeed  on  the 
total  systemic  dose  for  various  load  distributions.  The  results  for  five 
load  distributions  are  shown  in  Table  1 and  Table  2 for  windspeeds  of 
6.87  cm  sec-1  and  100  cm  sec"1  respectively. 

Also,  using  the  same  five  initial  load  distributions,  the  pro- 
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gram  was  applied  for  5,  10  and  20  slices  in  each  of  two  layers  and  times 
of  1.0  and  0.1  minutes  combined  as  shown  in  Table  3.  The  systemic  dose 
and  percentage  of  initial  mass  are  also  recorded  in  Table  3 to  show  the 
effect  of  varying  the  number  of  clothing  slices  and  the  time  increment 
size  on  the  results.  In  these  tests,  the  decomposition  was  set  to  zero 
so  that  the  total  mass  of  agent  should  theoretically  remain  constant. 

The  mass  deficiency  indicated  in  Table  3 is  then  only  due  to  approxi- 
mations caused  by  discretizing  the  skin  clothing  system. 

7.  DISCUSSION 

The  results  shown  in  Table  1 and  Table  2 indicate  that  the  wind- 
speed  over- the  clothing  surface  has  a considerable  effect  on  the  systemic 
dose  for  all  five  loading  configurations  tested.  The  largest  effect  occurs 
when  the  initial  load  is  placed  at  the  surface,  and  the  effect  decreases 
as  the  initial  load  is  distributed  closer  to  the  skin.  Also,  the  systemic 
dose  is  increased  as  the  initial  load  is  distributed  closer  to  the  skin, 
as  one  would  intuitively  expect. 

The  results  in  Table  3 for  various  time  step  sizes  and  number  of 
clothing  slices  indicate  that,  for  practical  purposes,  even  the  most  coarse 
discretization  produces  sufficient  accuracy.  Another  method  of  discreti- 
zation was  tried  for  the  slices  at  liquid-vapour  interfaces.  Instead  of 
assuming  that  the  diffusive  resistance  started  exactly  at  the  interface  as 
was  done  in  the  present  model,  except  for  the  skin-clothing  interface,  the 
assumed  resistance  included  half  that  of  the  liquid  containing  slice.  The 
difference  in  results  between  these  two  methods  was  found  to  be  less  than 
one  percent,  for  the  test  data  used  to  obtain  the  results  in  Table  1. 
Therefore,  either  method  produces  sufficiently  accurate  results  for  prac- 
tical purposes. 

8.  CONCLUSIONS 

A mathematical  model  and  an  associated  computer  program  have  been 
produced  to  calculate  the  systemic  dose  for  an  initial  liquid  loading  on  or 
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in  the  clothing  of  a clothed  skin  system.  The  program  can  be  applied 
to  practical  problems  for  any  liquid  agent,  when  the  diffu  ve  and 
absorptive  properties  of  a given  system  are  known. 
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TABLE  1 


COMPARISON  OF  SYSTEMIC  DOSES  FOR  VARIOUS 
DISTRIBUTIONS  OF  1.2  mg  cm"2  LOADING  AND 
WINDSPEED  OF  6.87  cm  sec-1 


INPUT  DATA 

RESULTS 

Loading  mg  cm-2 

Decay 

MS 

... 

DELT 

Systemic  Dose 

Surface 

Top 

Layer 

Bottom 

Layer 

min"1 

No. 

min 

at  29  Hours  mg  cm"2 

1.2 

0 

0 

0.003 

20 

0.1 

0.01297 

0.8 

0.4 

0 

0.003 

20 

0.1 

0.01359 

0.4 

0.8 

0 

0.003 

20 

0.1 

0.01417 

0 

1.2 

0 

0.003 

20 

0.1 

0.01473 

0.4 

0.4 

0.4 

0.003 

20 

0.1 

0.01541 
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TABLE  2 


COMPARISON  OF  SYSTEMIC  DOSES  FOR  VARIOUS 
DISTRIBUTIONS  OF  1.2  mg  cm-2  LOADING  AND 
WINDSPEED  OF  TOO  cm  sec"1 


INPUT  DATA 

RESULTS 

Loading  mg  cm-2 

Decay 

MS 

DELT 

Systemic  Dose 

j Surface 

1 

Top 

Layer 

Bottom 

Layer 

min-1 

No. 

min 

at  29  Hours  mg  cm-2 

! 

i 1-2 

0 

0 

0.003 

20 
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